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ABSTRACT 

We develop the theory of interstellar scintillation as caused by an irregular plasma 
having a power-law spatial density spectrum with a spectral exponent of P = 4 corre- 
sponding to a medium with abrupt changes in its density. An "outer scale" is included in 
the model representing the typical scale over which the density of the medium remains 
uniform. Such a spectrum could be used to model plasma shock fronts in supernova 
remnants or other plasma discontinuities. We investigate and develop equations for the 
decorrelation bandwidth of diffractive scintillations and the refractive scintillation index 
and compare our results with pulsar measurements. We consider both a medium con- 
centrated in a thin layer and an extended irregular medium. We conclude that the (3 = 4 
model gives satisfactory agreement for many diffractive measurements, in particular the 
VLBI meaurements of the structure function exponent between 5/3 and 2. However, 
it gives less satisfactory agreement for the refractive scintillation index than does the 
Kolmogorov turbulence spectrum. The comparison suggests that the medium consists 
of a pervasive background distribution of turbulence embedded with randomly placed 
discrete plasma structures such as shocks or HII regions. This can be modeled by a com- 
posite spectrum following the Kolmogorov form at high wavenumbers and steepening at 
lower wavenumbers corresponding to the typical (inverse) size of the discrete structures. 
Such a model can also explain the extreme scattering events. However, lines of sight 
through the enhanced scattering prevalent at low galactic latitudes are accurately de- 
scribed by the Kolmogorov spectrum in an extended medium and do not appear to have 
a similar low-wavenumber steepening. 

Subject headings: radiative transfer — ISM: general — pulsars: general 
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1. INTRODUCTION 

It is well known that radio waves propagating in 
the intcrstchar medium (ISM) arc scattered by the ir- 
regularities in the Galactic electron density. The scat- 
tering in turn gives rise to a number of observable phe- 
nomena. Among others, these include angular broad- 
ening and intensity fluctuations (both in time and fre- 
quency) of compact radio sources. While a nuisance 
in many radio astronomical observations, these phe- 
nomena can be used to investigate the nature of the 
irregularities in the interstellar plasma density. These 
density irregularities are in turn believed to follow 
the fluctuations in the interstellar kinetic and mag- 
netic energies. Ideally one would like to invert obser- 
vations of radio scintillation and scattering to deter- 
mine the statistics of the plasma density. As noted by 
Narayan (1992) this inverse problem is not well posed 
and one must rely on modeling methods. A com- 
plete prediction of scintillation observables requires 
an a priori knowledge of both the form of the spatial 
power spectrum of the electron density fluctuations 
and the distribution of the scattering material along 
the line of sight. For a given profile of the distribution 
of the scattering material, one may compare observa- 
tions and predictions and so constrain the functional 
form of the spectrum. The power spectrum provides 
useful insight into the physics of the plasma irregu- 
larities. Hence a knowledge of the form of the density 
spectrum becomes central in both predicting scintil- 
lation phenomena and understanding the physics of 
the interstellar plasma. In this paper, we revisit the 
investigation of the form of the density spectrum. 

A commonly used model for the density spectrum 
has been based on a power-law model with a large 
range between "inner" and "outer" scales (e.g. Rickett 
1977): 

Here q is the magnitude of the three-dimensional 
wavenumber q. 

(z) denotes the strength of fluctuations (with a 
weak dependence on distance z). /3 is the spectral ex- 
ponent, and K^^ = Li and = are the inner and 
outer scales respectively. For Ko « q « Ki, we ob- 
tain the sim,ple power-law model: Pnc{q) = C%^Q~^- 
Armstrong, Rickett, & Spangler (1995) have con- 
structed an empirical density spectrum by combin- 
ing radio scintillation observations in the local ISM 



(< 1 kpc) with measurements of the differential Fara- 
day rotation angles and large-scale electron density 
gradients. They have shown the power spectrum to 
be consistent with a simple Kolmogorov power-law 
model (/? = 11/3) over an astronomical 10 orders of 
magnitude in wavenumber scale (10~^^m~^ < q < 
10~^m~^). The Kolmogorov spectrum in density sug- 
gests a turbulent cascade in the magnetic and kinetic 
energies. This has lead to several theoretical investi- 
gations of the generation and maintenance of hydro- 
magnetic turbulence in the ISM (e.g. Pouquet, 1978; 
Higdon 1984 & 1986; Biskamp, 1993: Sridhar & Gol- 
dreich 1994; Goldreich & Sridhar 1995 and 1997). 
However, the Armstrong et al. study combined ob- 
servations from many lines of sight and the scatter 
among them leaves a substantial uncertainty in the 
exponent /?. A list of symbols is given in Table 3. 

In spite of the positive evidence for the simple Kol- 
mogorov spectrum, substantial observational incon- 
sistencies remain. For instance, long-term refractive 
intensity scintillations of some pulsars have modula- 
tion indices as much as a factor of 2 larger than pre- 
dicted by the simple Kolmogorov model (cf. Gupta, 
Rickett, & Coles 1993). Other discrepancies are re- 
vealed in the diffractive dynamic spectra of pulsars. 
On some occasions, periodic fringes are observed, 
which are not predicted by the simple Kolmogorov 
model (cf. Roberts & Abies 1982; Cordes & Wol- 
szczan 1986; Rickett, Lyne, & Gupta 1997; Gupta, 
Bhat & Rao, 1999); in addition, sloping bands in the 
dynamic spectra often persist longer than predicted 
by the model (cf. Gupta, Rickett, & Lyne 1994; Bhat, 
Rao & Gupta, 1999b Bhat, Gupta & Rao, 1999c). 
Further for some pulsars, the decorrelation bandwidth 
has larger amplitude variations than predicted for the 
Kolmogorov spectrum (Bhat, Gupta & Rao, 1999c). 
Such anomalies suggest the presence of large refrac- 
tive structures giving rise to the focusing and defo- 
cusing of the scattered ray bundles. The interference 
of the ray bundles can also explain the occasional 
fringes observed in the dynamic spectra of some pul- 
sars. Given the relatively frequent occurrence of such 
events, one can ask whether they can be considered as 
mere occasional anomalies or should be considered as 
a widespread phenomenon intrinsic to the spectrum 
on a grand scale. 

The observational inconsistencies suggest the need 
for an enhancement in the power on the large "re- 
fractive" (10"'^"'^ m to 10^^ m) spatial scales relative to 
the power on the small "diffractive" (10^ m to 10* m) 
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scales. There are several means by which this ratio 
may be enhanced. One is to include the inner scale 
cut-off in the density spectrum, which reduces the 
power at small scales (Coles ct al. 1987); these authors 
proposed inner scales of 10^ — 10^ m., though this does 
not correspond to any obvious physical scale. Physi- 
cally, the inner scale corresponds to the scale at which 
the turbulent cascade dissipates and becomes a source 
of heating for the plasma (Spangler, 1991). The value 
of the inner scale is largely unknown. Using different 
methods, Spangler & Gwinn (1990), Kaspi & Stine- 
bring (1992), and Gupta et al. (1993) have reported 
values for the inner scale ranging from 10'* to 10^ 
meters. In proposing the smaller values, Spangler & 
Gwinn (1990) argued that the inner scale is the larger 
of the ion inertial length, Li = Va/^i (where Va is the 
Alfven speed, and f2j is the ion cyclotron frequency), 
and the ion Larmor radius, = Vth/^i (where vth 
is the ion thermal speed); they obtained parameters 
for the warm ionized medium in reasonable agreement 
with observations. In a recent discussion, Minter & 
Spangler (1997) have suggested ion- neutral coUisional 
damping and wave-packet steepening as possible dis- 
sipation mechanisms for the turbulence in the diffuse 
ionized gas, which would make the mean- free path 
for ion-neutral collisions a possible value for the in- 
ner scale. However, this is thought to be larger than 
the maximum values proposed to explain the obser- 
vations, making it a less convincing dissipation mech- 
anism. Observationally, the Kolmogorov model with 
a large inner scale predicts refractive modulation in- 
dices consistent with pulsar measurements (Gupta et 
al. 1993). It has also been proposed to explain the 
occasional periodic fringes, with an inner scale on the 
order of the Fresnel scale (Cordes, Pidwerbetsky, & 
Lovelace 1986; Goodman et al. 1987). However, Rick- 
ett et al. (1997) reported a fringe event for the pulsar 
B0834-I-06 that could not be explained as the effect of 
a large inner-scale spectrum. The event requires sim- 
ilar conditions to those needed to explain the extreme 
scattering events (Fiedler et al. 1987). 

Another way to enhance the ratio of the power 

between the refractive and diffractive scales in the 
spectrum is to steepen the spectrum — with spectral 
exponents /3 > 4 (Blandford & Narayan 1985; Good- 
man & Narayan 1985; Romani, Narayan, & Blandford 
1986). While power-law spectra with f3 ~ 11/3 have 
a turbulence connotation, spectra with 4 < /3 < 6 
might involve some forms of turbulence, but are also 
consistent with a distribution of non-turbulent struc- 



tures with a range of spatial scales. Such steep spec- 
tra with 4 < /3 < 6 predict refractive modulation 
indices close to unity (Goodman & Narayan 1985), 
which is substantially larger than the range 30 to 40 
% observed from the nearby pulsars. On this basis 
Rickett & Lyne (1990) and Armstrong et al. (1995) 
have rejected spectra steeper than 4 for the inter- 
stellar plasma. However, the special case of /? = 4 
has been given little attention. We can conceive of a 
power-law model with spectral exponent /3 = 4 given 
by: 



Pn, {q; z) 



(2) 



Hereafter, we refer to this model as the "/? = 4 
model." Blandford & Narayan (1985) briefly dis- 
cussed this special case without including a cut-off at 
low wavenumbers. It is interesting to note that even 
though the Kolmogorov spectral exponent /3 = 11/3 
is very close to 4, the /3 = 4 model has a very differ- 
ent physical implication. Its physical origin has been 
rarely discussed, and it has not been formally com- 
pared with observations. Physically, this spectrum 
suggests the random distribution in location and ori- 
entation of discrete discontinuous objects across the 
line of sight. An "outer" scale, Lo = k^^, is in- 
cluded to account for the typical size of such objects. 
The "inner" scale here would correspond to the scale 
of the sharpness of a typical discontinuity. We as- 
sume the "inner" scale to be smaller than the diffrac- 
tive scale of scintillations, and hence has no signifi- 
cant effect on the scintillations. The /? = 4 model 
could characterize stellar wind boundaries, supernova 
shock fronts, sharp boundaries of HII regions at the 
Stromgren radius, or any plasma "cloud" with suffi- 
ciently sharp boundaries (transition regions shorter 
than the diffractive scintillation scale) which may 
cross the line of sight. Note that, though turbulence 
is not necessarily implied by the model, strong turbu- 
lence which has steepened to form shocks would also 
be described by the /? = 4 model. In these models a 
single discrete object crossing the line of sight could 
also explain the "extreme scattering events" observed 
in the flux density variations of extra-galactic sources 
(e.g. Fiedler et al. 1987). We note that the analysis of 
such events have been primarily based on geometri- 
cal optics involving a single "cloud" (Goodman et al. 
1987; Cordes ct al. 1986; Cordes & Wolszczan 1986; 
Roberts & Abies 1982; Ewing et al. 1970). Our anal- 
ysis of the (3 = 4: model includes the "wave optics" ef- 
fects when the line of sight passes through very many 
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such clouds. The ISM is assumed to consist of a ran- 
dom assembly of discrete structures with abrupt den- 
sity steps which may be independent of each other. If 
the /9 = 4 model were compatible with all of the scin- 
tillation observations, it would remove the implica- 
tion of interstellar plasma turbulence with an inertial 
range spanning as much as 10 orders of magnitude 
in scale, which has become the canonical model for 
ISS phenomena. It would, however, still be consis- 
tent with turbulence that has steepened into shocks. 

We start with a simple derivation of the (3 = 4 
model in section ^ In section |^, we give equations for 
the wave structure function for the /? = 4 model and 
compare predictions In section ^, the variation of the 
diffractive decorrelation bandwidth with frequency is 
used to test the (3 = 4 model against the simple Kol- 
mogorov model and the Kolmogorov model with an 
inner-scale (hereafter, the "inner-scale" model). In 
section ||, the observed variation of the refractive scin- 
tillation index with the normalized diffractive decor- 
relation bandwidth is compared with theoretical pre- 
dictions for the simple Kolmogorov, inner-scale, and 
(3 = A models. We use theoretical results from a previ- 
ous paper (Lambert & Rickett 1999; hereafter, LR99), 
in which we developed the theory of diffractive scin- 
tillations in a medium modeled by these spectra, and 
we add details of the theory for refractive scintillation 
in Appendix A. In Appendix B, we discuss the use of 
the square of the second moment to approximate the 
intensity correlation function for diffractive scintilla- 
tion with /3 = 4. In section ^ we give a discussion and 
our conclusions. 

2. THE SPECTRUM FOR RANDOM 
DISCONTINUITIES IN DENSITY 

In this section, we give a simple derivation of the 
/? = 4 model. Consider first a random distribution 
in space of identical plasma structures (blobs). The 
resulting electron density may be described as the 
convolution of a three-dimensional Poisson point pro- 
cess with the density profile of one individual blob. 
This is the spatial analog of shot noise, in which each 
charge carrier contributes the same temporal profile 
of current arriving at randomly distributed times. In 
the wavenumber domain, the power spectrum of the 
electron density is then given by the spectrum of the 
Poisson point process multiplied by the squared mag- 
nitude of the Fourier transform of one blob. The spec- 
trum of the Poisson point process is a constant equal 



to the number density, no, of the blobs in space (Pa- 
poulis 1991). Thus the power spectrum of the medium 
follows the same shape as that of an individual blob. 
A similar description applies to water droplets in a fog 
and many naturally occurring media, as discussed by 
Ratcliffe (1956). If the blobs are asymmetrical and 
randomly oriented, one must also average the spec- 
trum over the possible angles of orientation. 

The basic feature of structures with discontinu- 
ous boundaries is that their power spectrum has an 
asymptotic behavior at large wavenumbers that varies 
as (wavenumber)"^. This is the spatial analog of 
the idea that the temporal spectrum of any pulse 
shape with an abrupt rise or fall has a high-frequency 
asymptote as (frequency)"^. The prime example is a 
rectangular pulse, for which the spectrum is a sinc- 
squared function, whose high-frequency envelope fol- 
lows this law. The simplest spatial example is spher- 
ical blobs of radius a with uniform density /o inside 
and zero outside; clearly, an added uniform back- 
ground density does not affect the power spectrum. 
The three-dimensional spatial Fourier transform of an 
isotropic function /(r) of radial distance r is given by 
(Tatarskii 1961): 



F{q) = 



1 



27r29 



/(r) sm{qr)r dr . 



(3) 



where q is the magnitude of the wavenumber. For any 
function /(r) that falls to zero beyond a radius a, the 
integral will vary as 1/g in the high- wavenumber limit 
qa 1, making F{q) oc 1/q^. For the spherical blobs 
of size a we have: 



F{q) 



fo ( a cos{qa) 



sm{qa) 1 



(4) 



For the squared magnitude of the Fourier transform, 
for very small and large values of qa: 



\F{q)\' 



/o^aV(367r^) , qa « 1 , 



a2/2/(8^4) , qa»l. 



(5) 



The spectrum of the electron density is then no|F(g)p. 
We can approximate this by equation (^, which ex- 
actly matches equation ^ at small and large values 
of qa, with Cjf = n^a^ /q / (Stt^) and = Lq = 
(2/9)0-25a. 

There are evident generalizations to make the model 
more like a real medium. The spheres could have a 
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probability distribution for their radii, which would 
weight the average of the square of equation (^) . The 
spheres could have smoothly varying density inside 
and an abrupt boundary at radius a. As with the 
pulse example, the high-wavenumber behavior would 
not change. More realistic models would include 
anisotropic structures with random orientations and 
a distribution of scales. For example, consider a cir- 
cular cylinder as a model for an anisotropic blob, 
with a particular radius and thickness. Its Fourier 
transform can be computed with respect to its ma- 
jor axes. This can then be transformed to tilted axes 
and the result averaged over an isotropic distribution 
of angles of tilt. Evidently the result will also be 
isotropic; we have solved this approximately and find 
a spectrum with the same high-wavenumber asymp- 
totic form as equation (||), where the "outer scale," 
I/O, is the smaller of the thickness and radius. This 
result could be averaged over a probability distribu- 
tion of radii and thicknesses and a similar conclusion 
reached, with as a weighted average of the small- 
est dimension. Planar sheets at random angles are 
modeled as very fiat cylinders and provide a crude 
representation of a distribution of shocks. Another 
simplified model would be ellipsoids of different sizes, 
eccentricities, and orientations, and again we expect 
a similar result. 

We conclude that the /3 = 4 spectrum is an ap- 
proximate description of the power spectrum for a 
medium with a random distribution of discontinuous 
structures, in the limit for wavenumbers greater than 
the reciprocal of the smallest dimension across the 
structures. We note that for wavenumbers at or be- 
low about Ko — L^^ the functional form depends on 
the detailed shapes, though it will be much less steep 
than Thus equation (j^ will not be generally 

applicable at the lowest wavenumbers. Similarly, in 
practice there is a finite scale over which the den- 
sity jump occurs (e.g. shock thickness). This can 
be modeled by convolving the idealized discontinu- 
ous structures with a suitable narrow, say, Gaussian 
function which thus provides a very high-wavenumber 
cut-off beyond which the spectrum falls rapidly to 
zero. Consequently, we consider a spectrum with 
an inverse fourth power law between inner and outer 
scales, though these do not have the same connota- 
tions as inner and outer scales in a turbulent cascade. 
We do not discuss the physical origin of the supposed 
discontinuous plasma structures beyond the ideas for- 
mulated in the Introduction, such as supernova shock 



fronts, stellar wind boundaries or sharp boundaries of 
HII regions at the Stromgren radius. There needs to 
be no physical coupling between the individual blob 
structures, and there is no turbulent cascade implied 
by the /3 = 4 spectrum. However, we note that shocks 
developing from the steepening of very strong turbu- 
lence could also be described by the model, and so if 
successful, the model would not rule out such strong 
turbulence, though in this case the /3 — 4 range does 
not correspond to an incrtial cascade as does the 11/3 
spectrum. 

In the following sections, we give the necessary 
scintillation theory needed to develop equations for 
the scintillation observables corresponding to the /3 = 
4 model, which we then compare with ISS observa- 
tions. We concentrate on measurements which are 
sensitive to the form of the density spectrum and are 
relatively insensitive to the distribution of scattering 
material along the line of sight. 

3. PHASE STRUCTURE FUNCTION 

Central to the description of the second moment 
and the coherence function for intensity is the phase 
structure function. The longitudinal gradient of the 
phase structure function is related to the electron den- 
sity spectrum through (cf. Coles et al. 1987): 

^ Jo 

X [1 — Jo(ks)] fidK . (6) 

Here c is the speed of light, ro is the classical electron 
radius, and v is the radio frequency. For a plane wave 
incident on a scattering medium, the line of sight in- 
tegral of D'^{s, z; v) gives the structure function of the 
geometric optics phase, also called the wave structure 
function, D^i^s). For the spherical wave geometry ap- 
plicable to a pulsar at distance L, the separation s 
at the observer's plane is projected back to the point 
of integration, giving Ds{s) = /p D'^{sz / L, z;v) dz. 
The integrals are evaluated explicitly in equations 
(13), (17), and (21) of LR99 for a uniform scatter- 
ing medium with the simple power-law, inner-scale, 
and [3 — A spectral models, respectively. In the ex- 
tended scattering geometry, the field coherence length 
at the observer, sq, is defined by: £'s('So) = 1- 
For a screen of thickness Az at distance Zp from 
the pulsar, the wave structure function at the ob- 
server becomes Ds{s) = D'^{szp/ L, Zp;h')Az, which 
has an explicit dependence on the screen location. In 
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LR99, we found it useful to define a coherence scale 
So,scr for the screen by its phase structure function 
-^^(so.scr) = D'^{so^scr,Zp;v)^z = 1, which is inde- 
pendent of the screen location. Note that in applying 
equation (|6|), we are assuming that the screen thick- 
ness, Az, is much greater than the outer scale Lo and 
small compared to total distance L. 

Here, we repeat the equation for the phase struc- 
ture function of the /3 = 4 model in a thin layer 
(screen) . Substituting the /? = 4 model for the density 
spectrum into equation integrating, and multiply- 
ing by the screen thickness, Az, we obtain: 



[1 - {KoS)Ki{KoS)] 



(7) 



where Ki is the first order modified Bessel function of 
the second kind, and SM ~ C% Az is the scattering 
measure. Equation (Q) may be approximated by the 
following logarithmic expression: 



'In 



1 



(8) 



which matches the full Bessel expression at both large 
and small values of the argument. The detailed shape, 
near where the structure function flattens, is gov- 
erned by the shape of the low-frequency turnover in 
our model spectrum, equation (||); however, as noted 
above, the turnover shape depends on details of the 
density profiles in the plasma blobs which are not con- 
strained in our density model. 

The logarithmic structure function for the /3 = 4 
model has far-reaching consequences. As the outer 
scale becomes larger, the structure function approaches 
a square-law. However, due to the presence of the log- 
arithm, this occurs only slowly; that is, only at the 
limit as s goes to zero does the /3 = 4 model struc- 
ture function become exactly square-law. A square- 
law structure function for the medium is a conve- 
nient mathematical model, but as noted in LR99, in 
a real medium the structure function must eventu- 
ally saturate at very large separations. Nevertheless, 
the square-law structure function has been widely 
used since it provides a valid approximation for scales 
smaller than the reciprocal of the wavenumber at 
which the spectrum is cut-off. The fact that the 
/? = 4 model structure function deviates from a pure 
square-law, even for very large values of the outer 
scale, makes this model interesting and an indepen- 
dent investigation worthwhile. 
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Fig. 1. — Local logarithmic slope, a(s), of the phase 
structure function (upper panel) and the phase struc- 
ture function (lower panel) for the f3 — A model. We 
have set sq — 1000 km for all curves, and the solid 
lines correspond to different values of the outer scale, 
Lq. In going from the thinnest to thickest solid lines, 
the outer scales are: IQ-*, 10^ 10^, 10^, and 10" 
km, respectively. The thin dashed line corresponds to 
the simple Kolmogorov model with logarithmic slope 
of 5/3, and the thin dotted line corresponds to the 
square-law structure function. The points in the up- 
per panel correspond to measured local logarithmic 
slopes of the phase structure function obtained from 
VLBI (Spangler & Gwinn 1990). 
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In Figure |j, we have plotted the structure function 
of the /9 = 4 model and its effective local logarithmic 
slope, a{s), versus the transverse spatial lag, s, for 
various values of the outer scale, Lo (see equation 20 
of LR99). To illustrate the shapes, we have chosen 
So,scr = 1000 km for all curves. The thin dashed lines 
correspond to the simple Kolmogorov model with log- 
arithmic slope of 5/3, and the thin dotted lines corre- 
spond to the square-law structure function. The plot 
shows how slowly the local slope approaches 2; even 
for s/Lo ^ 10"^, the slope is 1.95-distinctly less than 
2. 

The points in the upper panel correspond to mea- 
sured local exponent of the wave structure function 
obtained from VLBI observations, as tabulated by 
Spangler & Gwinn (1990). In principle, this method 
of comparison provides a good test for the form of the 
density spectrum. These authors presented a similar 
plot corresponding to the inner-scale model, for which 
a{s) is close to 2.0 for s less than the inner scale 
and close to 5/3 for s greater than the inner scale, 
with a transition occurring over about a decade in 
s. They derived estimates of the inner scale ranging 
from 50,000 to 200,000 meters, for the highly scat- 
tered sources that they studied. Whereas there is 
some uncertainty in some of the estimates of a, there 
are several examples in which values greater than 5 /3 
are reliably observed (e.g. Trotter et al., 1998). As 
can be seen, there is also reasonable agreement with 
the /9 = 4 model. Assuming the model to be correct, 
a separate model-fitting to each observation yields an 
outer scale in the range Lq ~ 3 x 10^ to 10^ meters. 
These values are much smaller than the parsec scales 
of supernova remnants and interstellar clouds. How- 
ever based on this comparison alone, we are unable to 
discriminate between the inner-scale and P — A mod- 
els since they both provide equally satisfactory agree- 
ment with the observations. Hence further indepen- 
dent tests are needed, and in the following sections, 
we present comparisons with two other scintillation 
observables which in the end argue against the /9 = 4 
model as a universal description of the ionized ISM. 

4. DIFFRACTIVE DECORRELATION 
BANDWIDTH 

The diffractive decorrelation bandwidth, Avd, is 
perhaps the easiest ISS observable to measure. It is 
the frequency-difference for a decorrelation to 50% of 
the correlation function for diffractive intensity scin- 



tillations. These are usually recorded as a dynamic 
spectrum centered on a given radio frequency, u, for 
a pulsar. Observations of Av^ for one pulsar at a 
wide range of radio frequencies have provided an im- 
portant test of power-law models for the interstel- 
lar density spectrum on that particular line of sight. 
When the diffractive scale is far from the inner and 
outer scales, one expects Avd oc z^2/3/(/3-2)^ which for 
the simple Kolmogorov model is i/'* *, in both screen 
and extended medium geometries. In this section we 
re-examine the published measurements and compare 
them with theory for the three different spectral mod- 
els. Since the scaling laws are all so steep, in Figures 
§ and H, we have plotted the observations logarith- 
mically as Ai'd/i^'^- Overplotted are lines giving the 
theoretical predictions for our three spectral models, 
with the simple Kolmogorov model as a line of slope 
0.4. Before giving our conclusions, we must discuss 
the error bars and various corrections to the data and 
also the method for deriving the theoretical predic- 
tions. 

Observations 

Decorrelation bandwidth data, collected together 
by Cordes, Weisberg, & Boriakoff (1985; hereafter, 
CWB), for pulsars PSR B0833-45 (Vela), PSR B0329- 
-f54, PSR B1642-03, PSR B1749-28, and PSR 
B1933+16 are shown in Figures || and ^. We have 
also included recently measured points by Johnston 
et al. (1998) for the Vela pulsar. These are {v, Avd) = 
(8.4GHz, 13.9MHz) and (13.7GHz, 58.2MHz). The 
data are shown as solid circles if derived from dy- 
namic spectra and as solid stars if derived from pulse 
broadening measurements. The latter become easier 
to estimate at low frequencies, where the resolution 
bandwidth required for the former becomes too nar- 
row for an adequate signal-to-noise ratio. The con- 
version of an estimate of a pulse broadening time to a 
decorrelation bandwidth relies on the uncertainty re- 
lation 27rAi/dTrf Ci. Unfortunately, the "constant" 
Ci takes on different values for different geometries 
and spectral models. For the (3 — A model, it is also 
very weakly dependent on frequency. In Table 1 of 
LR99, we gave numerical values for Ci for the three 
spectral models with spherical waves in both a screen 
and an extended medium geometry. For a given geom- 
etry, the smallest Ci is for the Kolmogorov spectrum 
and the largest (by about 50%) is for the square-law 
structure function. Thus in Figures ^ and |^ we plot, 
as solid and open stars, the values converted using 
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Fig. 2. — Diffractive Decorrelation bandwidth versus 
frequency for PSR B0833-45 (Vela). For visual aid, 
a line of slope 4 has been subtracted from the log/log 
plots so that the frequency scaling for the square-law 
structure function model, indicated by the thick dot- 
ted line, would coincide with the abscissa. The filled 
stars correspond to decorrelation bandwidths com- 
puted from scattering broadening times using the Kol- 
mogorov value for the constant Ci , whereas the open 
stars correspond to the Ci from square-law structure 
function model. All lines are computed for the ex- 
tended medium, but would be nearly the same for a 
screen geometry. The thick dashed line corresponds 
to the simple Kolmogorov model. For the /3 = 4 
model (a), 4 solid lines have been drawn correspond- 
ing to parameter KqSo.ighz equal to 10~^^ (line clos- 
est to the square-law structure function model), 10~*, 
10^'^, and 10^^ (line flaring furthest from the simple 
Kolmogorov curve). For the inner-scale model (b), 5 
solid lines have been drawn for parameter KiSo,iGHz, 
ranging from 10^^ (line closest to the square-law 
structure function model) to 10^'^ (line closest to 
the simple Kolmogorov model), in equal logarithmic 
steps. 



Ci for the Kolmogorov and square-law structure func- 
tion models, respectively. We display the predictions 
for extended medium and note that the screen values 
are only about 20% less for each spectrum. Thus in 
comparing data with the theory for the /3 = 4 and 
inner-scale models, the data should lie between the 
extremes of the open and solid stars, depending on 
the outer and inner scales, respectively. Note that 
the value of Ci in equation (7) of the Taylor et al. 
(1993) catalog is more than 50% greater than values 
computed by LR99, since it concerns the mean pulse 
delay rather than the 1/e decay time of the pulse. 

The second correction to the data concerns the ef- 
fects of refractive scintillation at the higher radio fre- 
quencies, which are closer to the transition to weak 
scintillation. Gupta et al. (1994) gave a heuristic 
theory of the effect and expressions for the bias to 
the diffractive Ai'd due to refractive shifts; LR99 dis- 
cussed the effect from the work of Codona et al. (1986; 
hereafter, CCFFH). In Figures ||and|3[ we have plot- 
ted as open circles Av^ corrected for this effect using 
equation (D6) of Gupta et al. (1994), who also esti- 
mated the variability in estimates of Ai^d due to the 
refractive modulation. We used this to estimate an 
error bar on the open circles, which is typically larger 
than the error bars quoted by the original observers, 
plotted on the solid circles. 

In summary, the open circles and their error bars 
provide the best direct estimates of Av^ and the in- 
direct estimates lie between the open and solid stars. 

Theory 

The theory of diffractive scintillations has been dis- 
cussed by many authors; we will use the results of 
section 6 in LR99 for a point source (spherical waves) 
either in an extended scattering medium or with a 
screen at Zp from the pulsar and Zo from the ob- 
server (with Zp + Zo = L). Following usual practice, 
we assume the strong scattering limit, in which the 
frequency decorrelation function for intensity is the 
squared magnitude of the diffractive component of the 
two-frequency second moment of the field. The valid- 
ity of this approximation is examined in Appendix B, 
where we find that it introduces a small error for the 
(3 — 4: model, which however is negligible compared 
to the errors in the observations. The decorrelation 
bandwidth, Avd, can be defined in terms of a normal- 
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ized bandwidth, 



i^Vd = — 



where u — 



^scattC 



(9) 



Here v is the (geometric) mean of the two frequen- 
cies, and the parameter u determines the strength of 
scattering; c is the speed of hght and Zscatt is the 
effective scattering distance. For a uniform scatter- 
ing medium, Zscatt = L, and for a screen geometry, 
L. So is the field coherence scale 
where the appropriate phase structure function equals 
unity; for the uniform scattering medium, it is defined 
at the observer, and for a screen, it is so.scr which is 
independent of the distances to the pulsar and screen, 
and leads to User as the appropriate strength of scat- 
tering parameter. 

The computations presented in LR99 include a tab- 
ulation of the normalized bandwidth Vd for the various 
models under discussion, and associated plots of the 
intensity decorrelation function itself. For the sim- 
ple Kolmogorov model, a constant value for Vd is ob- 
tained regardless of the frequency or distance (0.773 
and 0.654 for the extended medium and screen geome- 
tries, respectively). Thus the frequency dependence 
of Ai/d reflects the frequency dependence of so , giving 
Avd oc :y2/3/(/3-2)^ CWB used this to estimate /3 from 
observations of Avd- We now compare the observa- 
tions with theoretical scaling laws for the two other 
spectral models. There are, however, some complica- 
tions. The quantity Vd is no longer exactly indepen- 
dent of frequency, and sq depends on both frequency 
and the outer scale, Lq, or the inner scale, Li, that 
parameterize the other two models. 

Consider the details for the discontinuity spectrum. 
In Figure ||, we plot Afd/i^Qn^ versus i^ghz as ob- 
served for the Vela pulsar. In Figure ^a, we have 
plotted theoretical curves for (3 — A model with vari- 
ous outer scales. We use equation (^ for the models 
and so need the variation of both sq and Vd with v. 

For cL screen, we determine so,scr 

from equation (^) 

for the (3 = 4 model, and eliminate the scattering 
measure SAI using the same equation for so,ighz at 
a reference frequency of 1 GHz. This gives a relation 
between So,scr/so,iGHz and i^ghz, with Lq/sq^ighz as 
a parameter. For an extended scattering medium, we 
use equation (21) of LR99 to determine sq at the ob- 
server, and the same method as for the screen to elim- 
inate SM. In Figure ^, we show the variation of Vd 
with Lo/so at a fixed frequency (computed by Lam- 
bert 1998 and described in LR99). This is combined 



with the so~'^GHz relation to obtain the relation of Vd 
to I'GHz- For the curves in Figure H, sq.ighz is deter- 
mined so that, for each value of Lg, the model Avd 
fits the measured values near IGHz. For most of the 
models of interest ^ sq, and then the logarithm 
functions vary much more slowly than the Sq term 
in equation (|8|). In which case there is nearly a lin- 
ear relation so/so,rof = v/vrcf- Furthermore, Figure^ 
shows that as so changes, the parameter Vd changes 
very slowly. Consequently, for a very wide range of 
the outer scale, a good approximation is Avd <x i/^, 
which reflects the fact that the underlying structure 
function then approaches a square-law behavior. 

Consider the plot for the Vela pulsar PSR B0833-45 
(Figure ^) for which the data are most extensive. Al- 
though we have computed curves and data corrections 
for both the screen and extended scattering geome- 
tries, we only show the extended medium results in 
Figure || because the two plots are so similar. In the 
comparison with the screen model, the star points 
are lowered about 0.08 vertical units (20%), and the 
/3 = 4 model curves are very slightly steeper functions 
of frequency. However, the scatter among the obser- 
vations is greater than the differences in the mod- 
els between screen and extended geometries. In the 
Vela plot, we see that the /3 = 4 model agrees some- 
what better with the observations than does the sim- 
ple Kolmogorov model. This pulsar is known to lie in 
a highly scattered region, and hence the presence of 
discontinuities, as incorporated in the /3 = 4 model, 
is perhaps reasonable. However, the conclusion is not 
strong, since the data are also reasonably fitted by an 
inner scale in the range 100 km to 500 km, as can be 
found from Figure 

For the other pulsars (Figure ||), the data show 
a stronger frequency dependence than predicted by 
both the simple Kolmogorov and /3 = 4 models. How- 
ever, the inconsistencies between the measurements at 
different frequencies are even greater, and better ob- 
servations are needed, before such comparisons could 
discriminate between the models. A recent series of 
measurements at 327 MHz (Bhat et al. 1999a, b,c) 
documents the variability of Avd, and a similar long 
sequence of such measurements at other frequencies 
is needed on the same pulsars, before reliable con- 
clusions can be reached from the frequency scaling 
observations. 
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Fig. 4. — Variation of the normalized difFractive decorrelation bandwidth, Vd- (a) The /3 = 4 model, upper curve 
for an extended medium geometry with parameter KqSq and lower curve for a screen geometry versus KoSo.scr- (b) 
Same plots for the inner-scale model. 



5. REFRACTIVE SCINTILLATION INDEX 

Slow variations in the flux of pulsars are caused 
by refractive interstellar scintillation (RISS), and can 
be characterized by the rms deviation in flux den- 
sity normalized by its mean, or scintillation index, 
mn. RISS is due to inhomogeneities in the interstel- 
lar electron density on scales much larger than those 
responsible for diffractive scintillation, as character- 
ized by the decorrelation bandwidth. There are now 
measurements of both phenomena on a substantial 
number of pulsars, from which we can constrain the 
density spectrum on each line of sight. The relation 
of the refractive scintillation index, rriR, to the nor- 
malized diffractive decorrelation bandwidth, Avd/i', 
depends on the ratio of the power in the density spec- 
trum at the large refractive scales (lO""^^ m to 10^^ m) 
to the power at the smaller diffractive scales (10^ m 
to 10** m). 

We compare pulsar measurements gathered from 
the literature with theoretical predictions by plotting 
mR versus Av^jv (cf. Rickett & Lyne 1990 and Gupta 
et al. 1993). Related tests have been made by Arm- 
strong et al. (1995), Bhat et al. (1999a) and Smirnova, 
Shishov, & Stinebring (1998). m^f measurements at 



610 MHz are listed in Table ^ and include the results 
of long-term monitoring observations by Smirnova et 
al. (1998), and older measurements near 100 MHz are 
listed in Table |[ The tables also include the decorre- 
lation bandwidth obtained from diffractive scintilla- 
tion, typically observed at a different radio frequency, 
scaled to the observing frequency for m/j. We used 
the Kolmogorov scaling law cx v'^ '^ , and we note 
that the minor differences that result from changes 
in the scaling law for other spectral models are in- 
significant in the comparison. The 100 MHz data 
consist of mfl measurements at 73.8 MHz, 81.5 MHz, 
and 156 MHz. The pulsars observed at 610 MHz are 
primarily located at distances >1 kpc, whereas the 
pulsars observed at 100 MHz are nearer at distances 
<1 kpc. Figures || and ^ show plots of the theoreti- 
cal and measured refractive scintillation index versus 
the normalized decorrelation bandwidth in separate 
panels for the two sets of measurements. Theoretical 
curves are also shown in the Figures for different spec- 
tral models in both the screen and extended medium 
geometries. The dashed line corresponds to the sim- 
ple Kolmogorov model [a = 5/3), and the solid lines 
in Figures ^ and ^ correspond to the /S — 4 and inner- 
scale models, respectively. 
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In section ^ we noted that variable refraction tends 
to decrease the measured decorrelation bandwidth 
and used an expression from Gupta et al. (1994) to 
apply a nominal correction to measured values. In 
the table of data from near 100 MHz, we have made 
use of the recently published measurements of Bhat et 
al. (1999a). They monitored the apparent decorrela- 
tion bandwidth near 327 MHz for 20 nearby pulsars, 
many of which have mp; measurements in Table ^. 
They found the bandwidth to vary by factors 3-5 and 
discussed the influence of varying refraction as a re- 
fractive bias. From their measurements, they derived 
a corrected decorrelation bandwidth, which we have 
scaled (using the v^-'^ scaling law) to the frequency 
at which mn was observed. These gave values typi- 
cally 2-5 times bigger than obtained from earlier mea- 
surements (e.g. Cordes 1986). The pulsar B0809-F74 
was not observed by Bhat et al. (1999a); however, 
the recent weak scintillation observations of this pul- 
sar by Rickett et al. (1999) similarly suggests that 
earlier decorrelation bandwidth measurements over- 
estimated the strength of scattering. The effect of 
these changes is to shift the plots to the right in Fig- 
ures ^ and H by about half a decade, compared to 
Figure 5 of Gupta et al. (1993). There is only one 
pulsar (B0329-I-54) in Table |l| common to the Bhat et 
al. (1999a) observations, and its decorrelation band- 
width has also been corrected for the refractive bias. 
Since the pulsars observed at 610 MHz were mostly 
more heavily scattered, the refractive bias correction 
is substantially smaller (see Gupta et al. 1994) and 
has been ignored. 

Theory 

The theory of refractive scintillations has been 
known since the 1970s and was applied to pulsar 
flux variations by Rickett et al. (1984). For exam- 
ple, Prokhorov et al. (1975) described how the mod- 
ulation index for intensity can exceed unity in strong 
scattering in their equations (4.53) et. seq. We use 
the notation of Coles et al. (1987) and confine our 
discussion to spherical wave sources, propagating in 
a scattering plasma which is either concentrated in 
a screen or extended uniformly between source and 
observer. The "low-frequency" approximation for the 
intensity covariance function is given by equation (10) 
of Coles et al. (1987), from which we obtain the nor- 
malized refractive variance, m^, by setting the spatial 



offset equal to zero: 

16^2 2^2^ r ^ , 

^ Jo Jo 
exp|-L^ D'^ [Krlj^h{x,y)] dyj x 

sin^ [O.bK^rp i^x{l — x)] ndndx, (10) 

where: 

h{x,y) = m:\n\y {I ~ x),x{l ~ y)] (11) 
rF,L = \JLc/2ttv. 



Here the transverse wavenumber is k = y '& ~^1y' 
and X = z/L, where z is the distance along the line 
of sight measured from the source. For the screen 
geometry, both the density spectrum and the gradient 
in the phase structure function D'^ are concentrated 
in a thin layer of thickness Az at distance Zp from the 
source and Zo from the observer. Thus the x and y 
integrations in equation ( p^ ) give x — y — Zp/L, and 
we find it useful to define an equivalent Fresnel scale 
rp.scr: 

^F.scr 

= yj Zcc/2'n:v where /L. (12) 

The K integration in equation (|o|) involves the prod- 
uct of the density spectrum, decreasing as a steep 
power of K, times the high-pass Fresnel filter oc 
up to fpscri times the low-pass exponential term 
that cuts off wavenumbers above s~^, where sr = 
UsctTy ,scr is the radius of the scattering disc. Hence 
in strong scattering [uscr > 1) we only consider the 

part of the Fresnel filter. Consequently, the inte- 
gration depends on the ratio Lq/ sr or Li/ sr. For an 
extended scattering medium the line-of-sight integra- 
tion softens the exponential cut-off in the k integra- 
tion, but the basic relationships remain the same. As 
others have noted the level of refractive scintillation is 
greater in an extended medium than in a screen with 
the same observed diffractive scintillation. 

The canonical spectral model in ISS studies of 
pulsars has been the simple Kolmogorov spectrum 
PnX^I^z) = C^Jz)g-"/3. By fitting this model to 
diffractive scintillation observations, observers have 
estimated the scattering measure, SM = C]^ {z)dz, 
toward many pulsars. When divided by the pulsar 
distance, this gives a line-of-sight average of Cj^ (z), 
which is found to vary greatly from one direction to 
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another and to increase dramatically for distant pul- 
sars seen at low galactic latitudes (e.g. CWB). Taylor 
& Cordes (1993) have developed a smoothed model 
for the Galactic plasma density distribution, which 
includes enhanced scattering in spiral arms and to- 
ward the Galactic center. However, there are also 
large random variations on much finer spatial scales, 
which would produce a scatter in a plot of mR against 
distance of dispersion measure. For a simple power- 
law spectrum both mn and ^v^jv depend on SM 
through a single strength of scattering parameter. 
Thus for the Kolmogorov spectrum the variation of 
mji with ^Vdjv is independent of distance or fre- 
quency and is given by a single dashed line in each 
plot in Figures |5| and |6|. If the scattering medium is 
uniform, the strength of scattering is u = t-p^l/sq. 
For the screen it becomes Ugcr — scr /so.scr, and in 
that case, the dashed theoretical line is independent of 
the location of the screen between the source and the 
observer. This point is substantiated in Appendix A, 
where the theory is laid out in more detail. According 
to this screen model, along each line of sight there is a 
scattering layer with a certain SM , which determines 
both rriR and Avd/v, but SM is not necessarily re- 
lated to the pulsar distance. The theoretical curve 
with SM as the variable is independent of where the 
layer is located along the line of sight. 

For the /3 = 4 and inner-scale models, the theo- 
retical curves for m/j versus ^v^jv depend on the 
extra parameter Lq or Li , respectively. Details of the 
theory are given in Appendix A, where the relevant 
parameters are shown to be Lo/so,scr and Li/so^scr- 
Since so,scr depends on frequency, SM, and distance, 
then mn also depends somewhat on frequency and 
distance. In order to fix the frequency dependence 
of the theoretical m^j values, we have separated the 
measurements into two groups (at 610 MHz and 
near 100 MHz). The distance dependence is dealt 
with by assuming that SM oc distance. Whereas 
this is clearly appropriate for the extended scatter- 
ing medium, it is less clear for the screen model, since 
the screen model supposes that SM is not necessarily 
related to distance. However, it is reasonable that, 
if on a long line of sight there is a single region that 
dominates the scattering, its SM value will statisti- 
cally increase with line-of-sight distance. Indeed the 
experimentally derived scattering measure increases 
faster with distance than if the medium were uniform 
(see Figure 1 of Cordes et al. (1991). 

In Figure we show the theoretical curves for mpj 



in the /3 = 4 model, which are relatively fiat where 
Lo > Sr. As u increases (Aiyd/i^ decreases), sr 
increases, and when sr > Lq the scintillations are 
suppressed. In a screen this occurs when l^Vd/v < 

(rp.scr/io)^- 

Inspection of Figure ^ shows that the measured niR 
values at 100 MHz are above the prediction of the sim- 
ple Kolmogorov model for the extended medium ge- 
ometry. These measurements can, however, be mod- 
eled with the /3 = 4 spectrum with suitable spe- 
cific choices of outer scale, Lo- The screen geometry 
is inconsistent with the measurements at 100 MHz 
for both the simple Kolmogorov and /3 = 4 models. 
For the 610 MHz data, the measured niR values are 
mostly in agreement with the prediction of the simple 
Kolmogorov model for the extended medium geome- 
try; however, they lie below the predictions of the 
/3 = 4 model. For the screen geometry, the measured 
niR values are above the simple Kolmogorov model 
curve. However, for the /? = 4 model, though less con- 
vincingly, an agreement with the measurements can 
be found by suitable choices of the outer scale, Lq. 
The most striking feature of the Figure is the good 
agreement with the uniform extended Kolmogorov 
model for most of the pulsars. This result, relies most 
heavily on the excellent data from Smirnova et al. 
(1998), who also note this agreement. We are per- 
suaded by this Figure that the (3 — A model cannot 
be viewed as a global alternative to the Kolmogorov 
spectrum. We also conclude that for the medium- 
distance pulsars measured at 610 MHz, niR does not 
agree with the screen geometry. It appears that even 
if the scattering medium is not uniform on scales of 
a few kpc, there is not a single region that dominates 
the scattering. Computations of niR due to a patchy 
distribution along the line of sight could test what dis- 
tribution would start to approximate the uniformly 
extended medium. 

Turning to a comparison with the alternative inner- 
scale spectrum in Figure ^, the 100 MHz observations 
could be explained by very large values of Li (lO^^m 
to lO^^m) in a screen configuration or by more mod- 
est, but still large, values of (10*m) in an extended 
medium. The latter was essentially the proposal made 
by Coles et al. (1987) for the nearby lines of sight. 
The 610 MHz points are relatively lower than those 
at 100 MHz. For an inner scale spectrum in a screen, 
large values of Li (lO^m to 10^'^m) would be neces- 
sary; and in an extended medium 60% of the points 
lie near the simple Kolmogorov model, with the rest 
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(a) 610 MHz, Screen 



(c) 100 MHz, Screen 





(h) 610 MHz, Extended Mediur 



(d) 100 MHz, Extended Mediur 



Logjo(Ai/yi/) 



Fig. 5. — Theoretical and measured refractive scintillation indices, m_R, versus the normalized diffractive decorre- 
lation bandwidth, /S.Vd/v. (a) and (b) show measurements at 610 MHz; (c) and (d) correspond to measurements 
near 100 MHz. The solid curves give the theoretical results for the (3 = 4 model with a range of outer scales from 
10^° m to 10^^ m, going from right to left in equal logarithmic steps. The dashed line corresponds to the simple 
Kolmogorov model in the strong scattering limit. Panels (a) and (c) give theory for a screen geometry; panels (b) 
and (d) for a uniform scattering medium. 
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Fig. 6. — Theoretical and measured refractive scintillation indices, m_R, versus the normalized diffractive decorre- 
lation bandwidth, ^v^jv. (a) and (b) show measurements at 610 MHz; (c) and (d) correspond to measurements 
near 100 MHz. The solid curves give the theoretical results for the inner-scale model with a range of inner scales 
from (corresponding to the simple Kolmogorov model) to The dashed line corresponds to the simple Kolmogorov 
model in the strong scattering limit. Panels (a) and (c) — screen geometry; panels (b) and (d) — uniform scattering 
medium. 
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requiring inner scale values (lO^m to lO^m). The 610 
MHz data in our analysis come from RISS observa- 
tions of 21 pulsars made by Stinebring and are de- 
scribed in more detail by Smirnova et al. (1998), who 
also compare the results with various models for the 
spectrum and spatial distribution of the electron den- 
sity. They note that 4 pulsars which are seen through 
known HIT regions or supernovae remnants show rela- 
tively elevated values for rriR, and find 3 other pulsars 
with similar behavior. These are the points which lie 
above the Kolmogorov line in Figure ^p. Their in- 
terpretation is that for these objects the scattering 
is concentrated into regions either near the pulsar or 
near the Earth, and that these regions are character- 
ized by an inner scale near 3 x lO^m. Their spectrum 
model is very similar to our inner-scale model, except 
that the cut-off is characterized by a steep power-law 
rather than by an exponential function. We find a 
somewhat larger numerical value for the inner scale 
required to match those objects. We also note that 
the theoretical m/j values are consistently higher for 
an extended medium, and so require a more modest 
inner scale. 

However, an important result of our analysis is 
an alternative explanation for the relatively high m/j 
values seen for some pulsars. We suggest that pul- 
sar lines of sight can pass through discrete clouds 
with increased plasma density on large scales, that 
steepen the low-wavenumber spectrum as opposed 
to cutting off the high-wavenumbers. A discontinu- 
ity spectrum (/? — 4) is one way that the spectrum 
can be steepened, but smoother structures on scales 
larger than lO^^m (~ 1 AU) would also boost the 
low-wavenumbers and could cause similar enhance- 
ments of refractive compared to diffractive scintilla- 
tions. Such enhancements are likely to be associated 
with HH regions or supernova remnants. A simi- 
lar idea was also discussed by Lestrade, Rickett, & 
Cognard (1998) in the context of extreme scattering 
events in pulsar timing measurements. To take this 
idea further, a composite between the Kolmogorov 
and the /3 = 4 or steeper spectra should be investi- 
gated, particularly one in which the line of sight is 
not uniformly weighted. Smirnova et al. (1998) also 
note that the strength of scattering increases very 
much faster with distance and dispersion measure 
than if the medium were statistically uniform. The 
distant pulsars in their sample are mostly observed 
at low galactic latitudes and so are subject to the en- 
hanced density and turbulence described by the Tay- 



lor & Cordes (1993) model. It appears that most dis- 
tant pulsars follow the uniform extended Kolmogorov 
model with inner scale smaller than lO^m. Thus al- 
though these lines of sight are subject to enhanced 
scattering in the inner galactic plane, the spectrum 
effectively follows the Kolmogorov law and the plasma 
is dispersed enough to approximate a uniform scatter- 
ing medium. There are presumably discrete refracting 
structures present along these lines of sight, but their 
contribution to the density spectrum is masked by the 
higher densities of the inner Galaxy, which still follow 
an apparently turbulent spectral form. 

6. SUMMARY AND CONCLUSION 

In this paper, the theory of the /? = 4 model for 
the electron density spectrum was derived for discon- 
tinuous density structures and compared with pulsar 
observations. A new feature of our analysis is the 
inclusion of an "outer scale" needed in any realistic 
model. The model is characterized by an effective ex- 
ponent a{s) of the structure function, which remains 
between 1.95 and 1.6 over a very wide range of Lq 
values (cf. Figure |l]). This at first seems a promis- 
ing explanation for the spread in the estimates of a 
derived from VLBI observations of the angular broad- 
ening profile, as observed on heavily scattered lines of 
sight. 

As discussed in section ^, from Figures || and ||, we 
find that the /3 = 4 model provides a somewhat bet- 
ter agreement with the measurements of the diffrac- 
tive decorrelation bandwidth versus frequency for pul- 
sar PSR B0833-45 (Vela) than does the simple Kol- 
mogorov model. This might arise from refractive scat- 
tering effects caused in the supernova remnant asso- 
ciated with the Vela pulsar. Four other pulsars with 
decorrelation bandwidths measured against frequency 
show an appreciably stronger frequency dependence 
than the predictions of both the simple Kolmogorov 
and (3 — 4, models. However, there are substantial in- 
consistencies among the measurements and better ob- 
servations are clearly needed, especially in view of the 
variability in Ai/^; documented by Bhat at al. (1999c). 

The predictions of the /3 = 4 model for the vari- 
ation of the refractive scintillation index with the 
diffractive decorrelation bandwidth are in partial agree- 
ment with the observations. As discussed in section ||, 
the values of mn measured near 100 MHz are above 
the prediction of the simple Kolmogorov model and 
had previously been explained as the effect of an inner 
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Table 1 

Refractive modulation index data at 610 MHz from the literature. 



Reference 

Rickett & Lyne (1990) 
Kaspi & Stinebring (1992) 



Stinebring et aL (1999) 
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B0136+57 


2.9 


610.0 


5.7x10-3 


0.15 


0.02 


B0329+54 


1.4 


610.0 


0.350 


0.37 


0.02 


B0525+21 


2.3 


610.0 


0.30 


0.31 


0.01 


B0531+21 


2.0 


610.0 


4.0x10^2 


0.32 


0.01 


B0736-40 


2.1 


610.0 


9.0x10"^' 


0.03 


0.01 


B0740-28 


1.9 


610.0 


3.5x10--^ 


0.13 


0.01 


B0818-13 


2.5 


610.0 


0.061 


0.23 


0.01 


B0833-45 


0.5 


610.0 


1.50x10-* 


0.24 


0.01 


B0835-41 


4.2 


610.0 


7.0x10""' 


0.21 


0.02 


B1641-45 


4.6 


610.0 


7.0x10-'^ 


~0.1 


0.05 


B1642-03 


0.5 


610.0 


0.770 


0.46 


0.04 


B1749-28 


1.5 


610.0 


6.0x10^2 


0.25 


0.02 


B1818-04 


1.6 


610.0 


4.0x10""' 


~0.1 


0.05 


B1859-03 


8.1 


610.0 


5.0x10"'' 


~0.05 


0.03 


B1911-04 


3.2 


610.0 


8.0x10""' 


0.22 


0.01 


B1933+16 


3.5 


610.0 


2.0x10-3 


0.24 


0.02 


B1946+35 


7.9 


610.0 


3.0x10-5 


~0.05 


0.03 


B2111+46 


5.0 


610.0 


3.0x10-3 


0.15 


0.01 


B2217+47 


2.5 


610.0 


0.20 


0.30 


0.01 



Table 2 

Refractive modulation index data at 100 MHz from the literature. 



Reference 


Pulsar 


Distance 
kpc 


1/ 
MHz 


Ai/ 
MHz 


rUr 




Cole et aL (1970) 


B0809+74 


0.31 


81.5 


1.38xl0"3 


0.45 


0.18 




B0834+06 


0.72 


81.5 


1.3X10-'"' 


0.39 


0.11 




B1919+21 


0.66 


81.5 


1.2xl0"3 


0.53 


0.18 


Helfand et aL (1977) 


B0329+54 


1.4 


156.0 


5.7x10-3 


0.35 


0.08 




B0823+26 


0.38 


156.0 


1.0x10-2 


0.46 


0.11 




B1133+16 


0.27 


156.0 


2.9x10-2 


0.50 


0.08 




B 1508+ 55 


1.93 


156.0 


7.0x10-3 


0.41 


0.07 




B1919+21 


0.66 


156.0 


1.4x10-2 


0.50 


0.07 




B2217+47 


2.31 


156.0 


6.45x10-* 


0.41 


0.11 


Gupta et aL (1993) 


B0329+54 


1.4 


73.8 


3.1x10-"' 


>0.15 






B0809+74 


0.31 


73.8 


8.90x10-"' 


0.34 


0.06 




B0834+06 


0.72 


73.8 


8.6x10-"' 


0.16 


0.03 




B0950+08 


0.12 


73.8 


0.17 


0.45 


0.05 




B1133+16 


0.27 


73.8 


1.6x10-3 


0.18 


0.02 




B1237+25 


0.56 


73.8 


1.8x10-3 


0.30 


0.06 




B1508+55 


1.93 


73.8 


3.8x10-4 


0.28 


0.09 




B1919+21 


0.66 


73.8 


7.5x10-* 


>0.21 
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scale, substantially larger than the values invoked by 
Spangler & Gwinn 1990. However, the measurements 
are also consistent with the /3 = 4 model with suit- 
able choices of the outer scale. For the 610 MHz data, 
most of the measured rrifi values are in good agree- 
ment with the prediction of the simple Kolmogorov 
model for the extended medium, and they lie be- 
low the curves of the /3 — 4 model. However, there 
are a significant number of good-quality observations, 
which lie somewhat above the Kolmogorov line. We 
suggest an alternative to a large inner scale on those 
lines of sight, that they pass through regions of en- 
hanced density, which cause enhanced refractive scat- 
tering; these regions must have less small-scale sub- 
structure than in a turbulent medium and could in- 
clude discontinuities. 

Based on the above considerations, we reject the 
/? = 4 model as a universal spectral model for the 
interstellar electron density fluctuations. The corol- 
lary is to strengthen the evidence for the Kolmogorov 
density spectrum, which in turn suggests a turbulent 
process in the interstellar plasma. However, the sim- 
ple Kolmogorov spectrum is not a universal model 
either, since it disagrees with several of the rrin ob- 
servations. Since the /? = 4 model provides reasonable 
agreement for many of these discrepant observations, 
we propose that enhancements in the large scale part 
of the spectrum (which need not be described by dis- 
continuities in density) occur on these lines of sight. 
With such enhancements causing the increase in re- 
fractive scintillation, there is no need to invoke the 
relatively large inner scales proposed by Coles et al. 
(1987). As proposed by Spangler & Gwinn (1990), a 
relatively small inner scale is then likely, controlled 
by the ion inertial length or Larmor radius. 

It appears that different spectral models need to 
be considered for different lines of sight. A widely 
distributed turbulent plasma with occasional large 
ionized structures that increase the effective average 
power density, Pn^ , at low wavenumbers (large scales: 
10"'^^ m to 10^'' m) is thus a model that needs further 
formal investigation. Very similar conclusions have 
been reached by Lestrade et al. (1998) and by Bhat 
et al. (1999b, c). This model could also explain the oc- 
casional "extreme scattering events" and episodes of 
fringes in dynamic spectra when a line of sight passes 
through a particular discrete density enhancement. 
New theoretical work is needed to quantify the ex- 
pected statistics of these propagation events. It is 
likely that numerical modeling will be necessary to 



model non-stationary scattering media. 
Acknowledgements: this work was supported by 
the NSF under grant AST-9414144. 



A. THEORY FOR REFRACTIVE SCINTIL- 
LATION INDEX 

Screen Geometry 

Here we derive expressions for the refractive mod- 
ulation index, mn, as a function of the normalized 
diffractive decorrelation bandwidth, Avd/i^, for the 
inner-scale and (3 — A models. The simple Kol- 
mogorov model is included as a special case where 
the inner scale goes to zero. A general integral ex- 
pression for the refractive scintillation index is given 
by equation (p^. With the medium concentrated into 
a thin layer (thickness Az at Zp from the source) and 
an assumption of isotropy in the density spectrum, 
the equation becomes: 



IGTr^r^c^A^ 



X sm 



/ 2 2 
2 / ''P.scr'* 



poo 

/ Pn, (k) exp [-D^{rl^^„K)] 
Jo 

(Al) 



K dn . 



For both the inner-scale and (3 ~ A models, the ex- 



ponential function in equation (Al) may be approxi- 
mated by the Gaussian (Coles et al. 1987): 



exp [-D^irl^^^^K)] ~ exp [-{tf. 



scr "-scr 



, (A2) 



where rp scr is the Fresnel scale and Ugcr is the strength 
of scattering, as defined in section ||. We now substi- 
tute models for the density spectrum P^ (k), given 
by equations (|l]) (with Kq = 0) and (^ and ob- 
tain relations between ttir and the scattering measure 
SM = C%^Az. To these we add the connections be- 
tween SM and so.scr, for each spectral model, from 
the equations for the phase structure functions (17), 
(11), and (21) of LR99. 

For the inner-scale model, we obtain: 



a2" r(l + a/2) 
s^,,„r(l-a/2) 

, Li 

1 + 



-^OjScr 



where 



q/(c 



(A3) 



(A4) 
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Table 3: Symbols Used 



Symbols Definition 



a, Q:{s) Exponent of separation s in structure function 

(3 Exponent of wavenumber in density power spectrum 

/3, /3' Vector spatial offsets in Appendix B 

c speed of light 

Ci Constant in uncertainty relation, 27rAudTd = C\ 

Cff Coefficient in electron-density wavenumber spectrum (m"^"/"^) 

Dff,,D'^ Wave/phase structure function and its longitudinal gradient Equation (6) 

Ds Wave/phase structure function for a point (spherical wave) source 

Az Screen or layer thickness 

Avd Diffractive Decorrelation Bandwidth 

e Fractional frequency difference \ki — fc2|/(fci + ^2) 

/('')i/o Radial profile function for electron density, interior density 

F{q) Three-dimesnsional Fourier Transform of /(r) Equation (3) 

g{T]),gi{r}) Intermediate functions oi t] — 

Jo Zero-order Bessel function of the first kind 

K Transverse (2-dimensional) wavenumber 

Ki, Ko High and low wavenumber cutoffs Equation (1) 

k = 27r/A, ki, k2 Radio wavenumbers 

km Geometric mean of fci , 

L Pulsar distance 

Li,Lo Inner and Outer scales 

A Radio wavelength 

iriR Scintillation index (rms/mean) for refractive scintillation Equation (10) 

mint Scintillation index due to scales intermediate between diffractive and refractive 

no Spatial number density of discontinuous objects 

u Radio frequency 

PnAq) Wavenumber spectrum for the electron density Equation (1) 

q Three-dimensional wavenumber 

Te Classical electron radius 

rp.L Fresnel scale ^^Ljk) Equation (11) 

rpjscr Fresnel scale for a screen {^J z^/k) Equation (12) 

s Transverse separation 

SM Scattering measure = Line of sight integral of C% 

So Field coherence scale 

sr Refractive scintillation scale 

Td Pulse decay time to 1/e 

u Strength of scattering parameter (= tf.l/so) Equation (9) 

Vd Normalized decorrelation bandwidth Equation (9) 

V4 Sum of wave structure functions Equation (Bl) 

X, y Normalized distances 

z, Zp Distance from pulsar along line of sight toward observer 

Zo Observer-screen distance 

Ze Effective screen observer-distance = ZoZp/ L 



Then substituting for SM we obtain an integral for 
m|j, which can be solved analytically using standard 
techniques (see e.g. Appendix 2 of Rickett 1973): 



(cosV')^"^^ cos(V'a/2) - 1 

[(Li/2so,scr)' + Ut 



4 ]"/2 
*scrj 



mR 



S0,s 



(2-a)/a 
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where Wgcr 



/so.scr, and ip is given by: 

2 



L, 



scr^sc 



(A6) 



We see that to/? depends on Uscr with £i/so,scr as 
a parameter. Similarly the difFractive decorrelation 
bandwidth, Ai/^/i', depends on Wgcr through equa- 
tion with as a parameter, which is a very slow 
function of £i/so.scr as discussed in section^. 

For the simple power- law model we set Li = 0, and 
in the strong scattering limit (uscr >> 1), equation 
(A6) gives 'ip ~ ''^scrj and we obtain the asymptotic 
expression: 

R 4 



('-f)^('-f) 



(A7) 



Putting this together for the simple Kolmogorov model 
with a = 5/3, we have the simple relation: mn = 
0.459(Aj/d/j^)°-^^^, plotted as the thick dashed lines 
in the screen plots in Figures || and ||. 

For the general inner-scale model, with Li non- 
zero, we show computed curves in Figure ^ We can 
recognize two regimes. Consider first the case of small 
Li. With sufficiently strong scattering, so,scr < Li < 



,scr } 



However, there is a 



complication in portraying this behavior in a plot ver- 
sus ts.Vdlv for fixed inner scale, because of the vari- 
ation of rp^scr with distance, which is not specified 
by the horizontal variable /S.Vd/v- We deal with this 
by obtaining an approximate scaling of rp.scr with 
Avd/v. When so,scr < -^i, equation (A3) relates 



-2 
'O.scr' 



' F,scr 



We now argue that 



on average SM cx pulsar distance cx rj, 
inating SM , we see that, r'F,scr — -scr 
scalings and so,scr < Li < rp.scr, we see: 

\2-a 



^ g^^, and ehm- 

0.5 



With these 



TO^ CX (Li/UscrrF,scr) " OC {L, / ,-,cfU^cr) 



1.5\2-a 



CX u J-^ OC (Az.rf/z/) 



0.25 



(A8) 



Here rp.rcf is the Fresnel scale for a "reference" pul- 
sar, and the last version in equation ( |A8D comes from 
using the Kolmogorov exponent a — 5/3, which gives 
the asymptotic slope of 1/8 for curves at the lower left 
of Fi gure Now with Li and rp.scr fixed, let 5o,scr 
increase (i.e. less scattering) until Li < so,scr < rp.scr- 
The inner scale is no longer important, and we get 
the same relation as for the simple Kolmogorov spec- 
trum: to|. cx Uscr*'^ cx {Aiyd/vY^^', visible where 
the curves steepen with increasing IS.Vdlv. 



For larger inner scales, the curves in Figure g show 
a pronounced peak. These occur for inner scales 
greater than the Fresnel scale and will be associated 
with focusing and caustics. Since our treatment only 
includes the first order term in the low wavenum- 
ber expansion, it is not reliable in the region of the 
peak. Goodman et al. (1987) have discussed caustics 
at length for the same inner-scale spectrum. They 
note that when so,scr < rp.scr < Li < sr, the scintilla- 
tion power spectrum starts to fill in at scales interme- 
diate between the diffractive and refractive wavenum- 
bers. Their equation (2.5.7) gives an estimate of the 
variance in this extra term as: 



TO;, 



2(^i/sfl)^ln(Li/so,scr) 

2(Li/Mscrrp,scr)^ln(ijUscr/rp,scr)-(A9) 



With fixed Li , mf^^^. increases much more steeply with 
decreasing Uscr than does to^. Thus as Uscrrp,scr 
decreases, there are fewer and fewer independent 
phase perturbations across the scattering disc. When 

"^scrrp scr ^ 

L^, focusing represented by to?„j ~ 
41n(Li/rp.scr) > 1 occurs. We note that the theo- 
retical mji versus Vd/^^ plots in Gupta et al. (1993) 
omit the focusing condition and are wrong for inner 
scales greater than the Fresnel scale. The effect of 
higher order terms in the low-wavenumber expansion 
has been studied by Dashen & Wang (1993). They 
obtain a more efficient expansion scheme that gives 
improved accuracy near the peak in scintillation in- 
dex. Nevertheless, it seems that a reliable prediction 
for the behavior near the peak in scintillations re- 
quires numerical evaluation. This becomes even more 
necessary in treating an extended scattering medium. 

We also note that the drop in mR as scattering 
gets weaker past the peak in Figure || is real. It repre- 
sents the fact that when a square-law structure func- 
tion applies for scales from so.scr up to the scatter- 
ing disc size, there is insufficient phase curvature and 
the scintillations remain weak even though Usa- > 1- 
In such circumstances the "scattering disc" is a mis- 
nomer, since an observer would see only a single an- 
gle of arrival, that could wander over a region of scale 

^scrrp.scr ■ 

Turning to the screen analysis of the (3 = 4 model, 
we relate so,scr to SM using equation (21) from LR99, 
where the phase structure function equals one. This 
gives the following analog of equation (A3) for the 
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inner-scale model: 



'O.scr 



In 



SO' 



(AlO) 

Substituting the ( 3 =^ 4 model for the density spec- 
trum in equation (|Al[) and letting rj ^ k^, we obtain: 



rrij 



1 



iv + i^lY 



exp[-(rF,scrWscr) ??] siu (rp ,,„77/2) dry (All) 

In order to evaluate this integral, we let g{rj) represent 
the integrand, and we let 51(77) be the same as g{ri) 
but with Ko set equal to zero. We can then rewrite 
the integral in equation (All) as: 



= 8{4,„.ln [l + (2L„/so,scr)']} ' X 
9i{l)dil~ I [giiri) - gmdA (A12) 

0^0 J 



The first integral can be evaluated analytically (see 
e.g. Gradshteyn & Ryzhnik 1965). In strong scatter- 
ing, the exponential term cuts off the oscillations of 
the sine term, which can be approximated by its ar- 
gument, and the 51(77) — 5(77) becomes negligible for 
values of 77 larger than k^. With these approximations 
we can also do the second integral. Putting these to- 
gether, we obtain: 



{ln[l + (2Lo/so,s 
tan-i(M-?)-2M; 



ln(l 



2C[(2 + C)e^Ei(C)-l]} 



(A13) 



where Ei is the exponential integral, C = (soKo)^uf(,j. — 
{srKo)'^- In the latter form sr — 7'F^scrUscr is the re- 
fractive scale (equal to scattering disk radius). Again, 
uiR is related to Aiyd/i^ through: Avd/v — Vd/u'^ci-- 
The resulting curves are shown in the screen panels of 
Figure |^. Consider the asymptotic behavior for 771,^ in 
equation ( [A13 ) as a function of Usa-- As the strength 
of scattering Ugcr increases, lS.Vd/v decreases (oc 



•'scr) 



In equation (All), the exponential term cuts off the 
integral (at l/s'^) before the oscillations of the sin^ 
Fresnel filter, which then approximates 77^rp If 
also C = {srKoY ^ li we can ignore in the denom- 
inator and the remaining 77"^ cancels the 77^ from the 
Fresnel filter, and the integral depends only on 1/ sr. 
In this approximation, 777,|j is then simply propor- 
tional to the slowly varying logarithmic term, which 



explains the relatively flat part of the curves in Fig- 
ure ^ under these conditions, in equation (A_13) the 
first two terms in the curly brackets sum to 2 and the 
last term is negligible. With fixed, now let User 
increase, making sj^ increase. Eventually C becomes 
greater than one when the scattering disk becomes 
greater than the outer scale. At this point the ex- 
ponential term cuts off the integral below Kg, where 
the spectrum fiattens. As itscr increases even further, 
the integral decreases steeply, causing the down-turn 
at very small IS.Vd/v- We again note that our ex- 
pressions rely on the first order of an expansion and 
will not be reliable near the peak in the scintillation 
index. However, there is not the same focusing con- 
dition that applied for very large inner scales. 

Extended Scattering Medium 

In order to obtain expressions for 771^^ for the two 
spectrum models in the uniform extended medium ge- 
ometry, we must complete the line-of-sight integrals in 
equation (^^ in addition to following the steps used 
in the screen geometry. For each distance x = z/ L m. 
the line of sight, there is also an integration over vari- 
able y in the exponential cut-off. If D'{s) oc s", this 
^-integration yields LD' [nr'^p ^x{\~'x)) / where 
rF,L = \J Lc/ [2t:v). This again provides a low-pass 
cut-off at the reciprocal of the radius of the effective 
scattering disk, where k ~ {usr-p,L)~^ ■ We define the 
scattering strength by us — r-p,L/sQg, with the field 
coherence scale sos, as in LR99, defined where the 
spherical wave structure function equals unity, mea- 
sured in the observing plane. For the other spectrum 
models there is not such a simple relation for the y- 
integration, but there is still an effective cutoff given 
by a similar equation. When the a:-integration is com- 
pleted, the effective Fresnel scale is actually smaller 
than rp.L due to averaging over the a/ x{1 — x). 

In analogy with the screen geometry, we make use 
of id entiti es similar to those given by equations (A3) 
and ( AlO ). For the extended medium, these identities 
are in turn derived from the wave structure functions 
for the inner-scale and P ~ 4 models (cf. LR99). The 
identities obtained thus will be similar to those for 
the screen geometry, except that for the inner-scale 
model, (1) there will be a factor of 3 on the right side 
of equation (A_3), and (2) the 1 in square brackets 
is replaced by [3/(1 + a)]"/^""^) ^ ^ g foj. ^j^g j^ol- 
mogorov exponent. For the /3 = 4 model, the only 
change will be a factor of 3 on the right side of equa- 
tion (AlO). We use all of the aforementioned idcnti- 
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ties for the inner-scale and 13 — 4 models and compute 
the X integral numerically, since it cannot be carried 
out analytically. 

The shapes of the curves bear a close relationship 
to the screen results, though the extended medium 
values generally lie above the associated screen values 
at the same Af^/t^. 



B. DifFractive Intensity Correlation Function 

The second moment of intensity is needed to de- 
scribe the fluctuations of intensity. Under strong scin- 
tillation conditions, separate forms can be used for 
refractive and diffractive fluctuations, since their spa- 
tial scales differ by several orders of magnitude. In 
LR99, as elsewhere in the ISS literature, the correla- 
tion of diffractive scintillations is approximated by the 
squared magnitude of the second moment of the field, 
leading to the simple result that the spatial scale of 
the diffractive scintillations is equal to the scale where 
the phase structure function equals unity (so). How- 
ever, Goodman and Narayan (1985) showed that for 
steep spectra {(3 > 4) this is no longer the case and 
the diffractive scale can be larger than sq. Here we 
examine this question for the (3 = 4: spectrum. We 
give the details for a phase screen with plane wave 
source, which are readily generalized to a spherical 
wave source. 

The two- frequency intensity cross-spectrum at wave- 
number K, for a screen at distance L is given by the 
Fourier-like integral equation (17) of CCFFH. This 
depends on the combination of structure functions V4, 
which for a plasma screen can be written as: 



Kn run 



where ki and /c2 are the two radio wavenumbers, km 
is their geometric mean, k is their arithmetic mean, 
ko = k'^Jk, € = \ki — k2\/2k, and is evaluated 
at km', is a spatial offset which is the variable of 
integration. 

Consider first the single frequency case fci = A:2 (e = 
0). In the limit of very large k, the first four struc- 
ture functions saturate and sum to zero. The last two 
are equal and V4 w 2Z?0 (/?'), which gives the simple 
diffractive limit mentioned above. This is the zero- 
order term of an expansion, which is obtained in terms 



of the sum of the first four terms as a small quan- 
tity. The zero order result requires full saturation, 
which requires nL/km ^ L^- In diffractive scintilla- 
tion K ^ 1/so; hence, the condition becomes that the 
refractive scale sn = L/{kmSa) > L^. Our concern 
here is to consider what happens when the diffractive 
K is not large enough for saturation of D^. For shal- 
low density spectra (/3 < 4) , small argument approxi- 
mations to the structure function follow an exponent 
(3 — 2 < 2, and the zero-order term gives a good ap- 
proximation even when sr < Lq. However, for steep 
spectra Goodman and Narayan (1985) showed that 
the leading term in the structure function follows a 
square law, which exactly cancels in V4; the result 
is that the high wavenumber limit depends on the 
next term in the structure function expansion, which 
yields a diffractive scale that is greater than the scale 
So (defined by the square law term). 

Now we consider the case for /3 = 4 model, when 
the scattering disc su is smaller than the outer scale 
Lo- Here we can approximate equation (^) by: 



Da 



In(sVsg) 



(B2) 



As for the steep spectra, the leading term in the struc- 
ture function follows a square law, which cancels when 
substituted into equation (Bl). V4 can then be ap- 



proximated for large k by expanding in (3'km/{KL). 
The result is: 



V4 



2P'^[\n{u'^) + 1 + 2cos^ 6 - \n{P'^/sl)] 



(B3) 



Here u is the strength of scattering defined in equa- 
tion (^; 9 is the angle between vectors /t and /3'. 
We note that for a Kolmogorov spectrum in the high- 
wavenumber limit, V4 also includes terms in cos^ 9, 
which would be accounted for in the higher order 
terms of the expansion. A result similar to equa- 



tion ( |B3| ) is given by Dashen & Wang (1993), though 
for a 1-dimensional phase screen. In considering the 
spectrum of intensity fiuctuations (eq. 17 of CCFFH), 
one can show that the dominant wavenumber is ap- 
proximately K ~ P'^^ With these substitutions in V4, 
which we then set = 2, we solve for (3'; this gives an 
approximate equation for the diffractive spatial scale 

Sd- 

2 _ 2 ln(4/s§^g) 
'° \n{u^) + 2 ■ 

Under the condition assumed in this approximation, 
SR ^ Lo, we find Sd > so. However, in practice the 
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ratio Sd/so never becomes large. With a large outer 
scale, say Lo = 3 pc, and typical observing conditions 
So ~ 10* m and u > 100, we find Sd ^ I.Tsq. Thus 
the diffractive scale could be 70% greater than sq and 
would slowly approach sq for smaller outer scales. 

Turning to the two-frequency intensity correlation 
(0 < e ^ 1) in the diffractive limit of large k, the 
results of CCFFH still apply. Namely, that the last 



two terms of equation (Bl) largely control the decor- 
relation versus frequency. They group the remain- 
ing terms into a filter that depends only on k and a 
smaller term that becomes the basis of the expansion. 
The filter term was discussed by LR99 and shown to 
be important only as the strength of scattering de- 
creases. It is the last two terms in V4 that determine 
the zero-order result, so we looked at the effect of 
the higher order terms. The quantity that we are 
ultimately concerned with is the cross-correlation of 
intensity at offset frequencies at the same observing 
point. This comes from the integral of the cross- 
spectrum. Equations (31) through (34) of CCFFH 
give the zero and first order terms of the cross spec- 
trum in terms of the spectrum of refractive index 
fluctuations in the layer. For the /? = 4 spectrum 
we reduced these to a sum of confluent hypergeo- 
metric functions, which can be explicitly computed. 
For a sample observing condition we found that the 
higher order terms for the cross-spectrum itself were 
significant compared to the zero-order term; however 
when integrated to give intensity cross-correlation, 
they only had a minor effect on the decorrelation 
bandwidth itself (<5% increase). The reason for this 
appears to be the dominant effect of the last two terms 
in the V4 summation with unequal frequencies. 

To summarize we find a modest (logarithmic) in- 
crease in the diffractive scale relative to the field co- 
herence scale So J but that this remains less than a 
factor of 1.7 for the likely ISS parameters. This is 
accompanied by a smaller increase in the decorrela- 
tion bandwidth relative to the calculations of LR99, 
which relied on the normal zero-order expansion at 
high wavenumbers. This small offset in the decorre- 
lation bandwidth is negligible compared to the mea- 
surement errors for the observations under considera- 
tion. We assume that the conclusions reached here for 
a screen would also apply for an extended scattering 
medium. 
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